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! Abstract 

i-£h ' At the order of 1/ mf , the B meson lifetimes are controlled by the hadronic 

o ■ matrix elements of some four-quark operators. The nonfactorizable mag- 

nitudes of these four-quark operator matrix elements are analyzed by 
QCD sum rules in the framework of heavy quark effective theory. The 
5-i ' vacuum saturation for color-singlet four-quark operators is justified at 

hadronic scale, and the nonfactorizable effect is at a few percent level. 
However for color-octet four-quark operators, the vacuum saturation is 
violated sizably that the nonfactorizable effect cannot be neglected for 
the B meson lifetimes. The implication to the extraction of some of the 
parameters from B decays is discussed. The B meson lifetime ratio is 
predicted as t(B~)/t(B°) = 1.09 ± 0.02. However, the experimental 
result of the lifetime ratio T(Ab)/T(B°) still cannot be explained. 
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Heavy hadron lifetimes provide us with testing ground to the Standard Model, 
especially to QCD in some aspects [1-3], because they can be systematically calculated 
within the framework of heavy quark expansion. Theoretically, if we do not assume the 
failure of the local duality assumption, the heavy hadron lifetime differences appear, 
at most, at the order of 1/mg [4]. Recent experimental results on the lifetime ratio of 
At baryon and B meson [5] showed some deviation from the theoretical expectation. 
This has drawn a lot theoretical attentions [6-11]. The current experimental values for 
the lifetime ratios which we are interested in are [5] 

t(B~) 

tW = L06±0 - 04 ' 

I|M = 0.79 ±0.06. 

This may imply that the O(-t-) contribution is not enough for the explanation of above 
heavy baryon and heavy meson lifetime difference. To the order of 1/mjj, the hadron 
lifetimes have been studied since mid-80s [11, 12, 4, 6, 7]. And the potential importance 
of the l/mf corrections has been pointed out. The parameterization of the hadronic 
matrix elements of four-quark operators which appear in the hadron lifetimes at the 
order of l/m\ is generally expressed as [7] 

{B\h^l - 7 5 )W7 M (1 - l5)b\B) = B 1 F^m% , 



(B\b(l- l5 )qq(l + l5 )b\B) = B 2 F 2 B m 2 



B > 



and 



<£|M1 - l^aqqY{l ~ l 5 )t a b\B) = e^ml , 
(B\b(l - lb )t a qq{l + lb )t a b\B) = e 2 F 2 B m 2 B , 

Flm B 



(2) 



■(AftlM 1 " 75)5^(1 -7 5 )6|Ai 



2m Ab " ^ 1 12 

j ■<A 6 |5(l- 7 6)??(l + 7fe)6|A*> = r>'"' 



(3) 



2m A / U| ^ — 24 ' 

where the parameters Bi, (i = 1,2), Fb, r and B should be calculated by some 
nonperturbative QCD method. In above equations, the renormalization scale is arbi- 
trary, and the parameters depend on it. It can be taken naturally at the low hadronic 
scale to apply heavy quark expansion. On the other hand, if the scale is taken at m b , 



parameter Fsirnb) is just the well-defined measurable physical quantity — B meson 
decay constant 



The QCD sum rule [13], which is regarded as a nonperturbative method rooted in 
QCD itself, has been used successfully to calculate the properties of various hadrons. 
In Ref. [8], in the framework of heavy quark effective theory (HQET), the baryonic 
parameters r and B have been calculated by QCD sum rule, r ~ 0.1 — 0.3, B ~ 1. 
For a complete analysis, the mesonic parameters Bi and e« should be also calculated 
from QCD sum rule. The four-quark operators, and hence B iy e«, r and B, are scale- 
dependent quantities when the QCD radiative corrections are included. It was proposed 
by Shifman and Voloshin [11] that at the low hadronic scale, the vacuum saturation 
approximation, namely Bi — 1 and 6j = 0, makes sense. However in this case, the 
measured lifetime ratio r(A b )/r(B°) cannot be explained [8]. There are some argu- 
ment, on the other hand, that the vacuum saturation maybe a poor approximation [9]. 
Especially from a naive large N c analysis, q's are about 1/N C ~ 0.3 [11, 7]. We will 
explore the violation of the vacuum saturation approximation in detail from QCD sum 
rules in the framework of HQET. 



Let us first consider the parameters B^. We construct the following three-point 
Green's function, 

= i 2 J ^e ifc ''^(0|T[^ , (4) 

where u — 2v ■ k, u' — 2v ■ k'; h$ is the b-quark field in the HQET with velocity v. 
And O denotes the color-singlet operators given in Eq. (2), 

O = bT iq qT 2 b , (5) 

with Ti = T 2 = 7 M (1 — 75) for B x and T x = 1 — 75, T 2 = 1 + 75 for B 2 . In terms of 
the hadronic expression, the parameter Bi appears in the ground state contribution of 

r>X), 

r>, a/) = fr (2A _J) (2 %_„ + resonances , (6) 
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where A = ms — mb. The resonance contribution will be simulated by the perturbative 
QCD contribution above some threshold energy due to the local duality assumption. 
On the other hand, this Green's function T°(u, uj') will be calculated in terms of quarks 
and gluons, that is to say, by the operator product expansion method of QCD. The 
essential feature of the QCD sum rule is that in the QCD calculation, the vacuum 
condensates of quarks and gluons have to be included. Practically the calculation is 
performed at 1 GeV scale or so, only a few terms of the condensates with lowest dimen- 
sions are important. Note that this calculation will be reinforced by the Borel trans- 
formation, and its consistency should be checked through the finding of the so-called 
sum rule window. This procedure will be explained in more detail in the calculation 
of parameters e^. 

The calculation of T°(uj,uj') in HQET is straightforward. The fixed point gauge 
[14] is adopted. The assumption that the four-quark condensates are factorizable is 
used. The violation of it will be discussed later. The dominant non- vanishing Feynman 
diagrams are shown in Fig. 1 where the double lines denote the heavy quark. However, 
all these diagrams are factorizable, namely they do not contribute any deviation from 
vacuum saturation. And this is true even when the 0(a s ) radiative corrections to 
these diagrams are included, simply because O is a color-singlet operator. Note that 
due to the same reason, there is no nonfactorizable gluon condensate contribution. 
The tadpole diagrams in which the light quark lines from the four-quark vertex are 
contracted have been subtracted. Generally, the nonfactorizable diagrams for a color- 
singlet four-quark operator are listed in Fig. 2. They would be the leading contribution 
diagrams which might have an 0(1) correction to B^ — 1. But it is easy to see that they 
are vanishing because of the special structure of Ti and T 2 in O (Eq. (5)). Hence, we 
see that the vacuum saturation approximation is valid at hadronic scale for color-singlet 
operators, Bi — 1, for i — 1, 2, through the leading order consideration. 

To what extent the nonfactorizable effect makes Bi deviate from unity? The next 
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possibility for non-factorization is to consider two gluon exchanges, like those in Fig. 
3. It is interesting to note that the four-gluon condensate < Oi 2 s G A > contribution 
to r°(co>,k/), like Fig. 3(b), is vanishing^. The perturbative two gluon exchange dia- 
gram, namely Fig. 3(a), can be the leading non- vanishing nonfactorizable contribution. 
(There are some other diagrams which have the same or less order of magnitude numer- 
ically, like Fig. 3(c).) This is a four-loop diagram, its contribution to Bi is estimated as 
^ « s (iGcV) ^2 which j s numerically around a few percent and is too small to be important 
for the hadron lifetimes. We obtain that 

^ = 1 + O(10- 2 ), (i = 1 , 2) . (7) 



The parameters e, 2), as we will see in the following analysis, have deviation 

from the expectation of the vacuum saturation at the hadronic scale which cannot be 
neglected for the B meson lifetimes. The procedure is similar to that for parameters 
B{. The three-point Green's function in this case is constructed to be 

r>,u/) = i 2 j dx^e^-^(0|^ , (8) 

where T is the color-octet operators given in Eq. (2), 

T = bY x t a qqY 2 t a b . (9) 

In the hadronic language, the parameter appears in the ground state contribution of 

r>x), 

T>, u') = ti {2K _* )(2 l_ U}l) + resonances . (10) 



In the calculation of T t (uj,uj'), all the condensates with dimensions lower than 6 

are retained. The dominant diagrams are found to be those given in Fig. 4, and 
2 This fact is gauge invariant and can be easily observed in the fixed point gauge. In this gauge, the 
light quark propagator with two condensate gluons attached vanishes if one end of the propagator is 
at space-time origin [14]. And in this gauge, the heavy quarks in Figs. 1-3 are free from interaction. 
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they are non- vanishing. The four-quark condensate (qq) 2 diagram vanishes. We have 
neglected the perturbative diagrams which are three- loop diagrams of order a s . That 
means the QCD radiative corrections to are not included. From general experience 
of QCD sum rule method, the condensate diagrams play more dominant role than 
the corresponding perturbative one. This neglect is expected viable. The condensates 
parameterize the non-perturbative effects which in the 1 GeV scale are still small 
enough to be treated as power corrections of 1/ou and 1/u/ in the operator product 
expansion. Because the perturbative contribution has been neglected, resonances in 
Eq. (10) will be also neglected due to the duality assumption. While the calculation 
can be justified if (— uj) and (-co') are large, however the hadron ground state property 
should be obtained at small (-co) and (— uj'). These contradictory requirements can 
be achieved by introducing double Borel transformation for uj and uj' . There are two 
Borel parameters T and T' corresponding to uj and uj', respectively. They appear 
symmetrically, so we take T = T' in the analysis. 



The sum rules for the parameters e$ are 

f e A ~ A/f 

( {* S GG) ~ \ f e^/ f (U) 

62 = -{-^ T+m ° {qq) )T^n^ B ' 

There is no gluon condensate contribution to e\. Numerically we use the following 
values of the condensates, 

(qq) ~ -(0.23 GeV) 3 , 

(a s GG) ~ 0.04 GeV 4 , (12) 

(gqa^G^q) = m 2 Q (qq) , m\ ~ 0.8 GeV 2 . 

For consistence, the HQET sum rule of the parameter Fg [15, 16] will be used. The 
result derived in Ref. [16] is 

= h r dm?e "" f - m - ■ <i3) 



where u c is twice of the continuum threshold, which was determined as u c ~ 2.0 ± 0.3 
GeV. The range of the Borel parameter T in Eq. (11) should be similar to that of 
Fb- This point can be obviously seen from the sum rules for Bi if it is written down 
explicitly. In that sum rules, if all the non-factorizable contributions are neglected, 
then the resulting E>i is unity only if the Borel parameter is the same as that of F B - 
Practically we take the window as 0.7 < T < 1.0 GeV. There is no A dependences for 
ej's, because they are canceled actually. The numerical sum rule results for e\ and e 2 
are given in Fig. 5. From the figures, we see that q's have no good stability in the 
window of T. This is because we have not included the perturbative diagram. Finally 
the results for e 1 and e 2 are obtained as 

Cl ~ -(4.1 ±2.2) x 10~ 2 , 

(14) 

e 2 ~ (6.1 ±3.5) x 10~ 2 . 

In spite of the large uncertainties in above equation, the numbers are significant for the 
B meson lifetime difference. We see that the vacuum saturation approximation for the 
color-octet matrix elements is indeed violated. The magnitudes of q can be as large as 
0.1. 

It is necessary to discuss the hypothesis on the four-quark condensate factorization. 
We have used it in above analysis for Bi and e«. While this is the working assumption 
for the usual QCD sum rule calculations, the violation of it may imply in our case, new 
contribution to the non-factorizable effect of the four-quark operator matrix elements. 
This contribution to B^ and q is estimated as 

4A/T 

6B i ~6e i = 6(qq) 2 -z-^- , (15) 
b B m B 

where S(qq) 2 denotes the deviation of the four-quark condensates from factorization. 
Arguments based on large N c expansion suggest that this approximation is good to 
within l/N 2 ~ 10% [17]. Numerically, considering the success of this assumption in 
the calculation of baryons by QCD sum rules [18], we take 5(qq} 2 to be less than 30% 
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of (qq) 2 . In this case, SBi and 8ei are smaller than about 10 2 which have no significant 
influence on our results for Bi and e,. 

It should be noted that the vacuum saturation and its violation we have analyzed 
above are at some hadronic scale, other than the scale of mj. Because we have been 
working in the framework of HQET, in which the natural scale is jihad <^ m b- The 
renormalization group evolution of the relevant operators was calculated in Refs. [19, 
7]. Information on parameters Bi and e« at scale m b is necessary to obtain the hadron 
lifetimes by using the analysis of Ref . [7] . There is a notation difference in this paper. In 
HQET, Fb given in Eq. (13) is in fact a scale-dependent quantity if the renormalization 
effect is considered. From this point of view, Fb in our previous analysis took value 
at hadronic scale, whereas that in Ref. [7] at scale of m b . Taking this difference into 
consideration, we have the following evolution relations by choosing a s (/ihad) — 0-5 
(corresponding to jihad ~ 0.67 GeV), 

Bi(m b ) ~ B^ad) - 0.24ei(/j,had) , 

(16) 

€i(m b ) ~ -0.05Bi(jj,had) + 0-72ei(nhad) ■ 

That is 

B x [m h ) ~ 1.01 ±0.01, B 2 (m b ) ~ 0.99 ± 0.01 , 

(17) 

ei(m 6 ) ~ -0.08 ±0.02, e 2 (m b ) ~ -0.01 ± 0.03 , 
from Eq. (14). 

To be more relevant, we discuss a more detailed parameterization for the four- 
quark operators proposed in Ref. [6] in the following. It is motivated for the model- 
independent determination of the hadronic matrix elements from the lepton spectrum 
in the endpoint region of semileptonic B decays. The parameterization is written as 

(Blhu-i 1 ~ 75)<?<hv(l - ls)b\B) = (v sing iv^v v - g S in g i9^)F^m 2 B , 

(1°) 

(-8|&7m(1 -ls)t a qqiv{^ -ls)t a b\B) = (voctV^Vv- g oct g„ u )F%m B , 
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where v S i ng i, g S i ng i and v oct , g oc t are the parameters. The tadpole diagram contributions 
of the matrix elements have been subtracted. The parameters are related to those 
of Eq. (2) as v sing i — Ag sing i = B\ and v oct — ^g oc t — £i- In the vacuum saturation 
approximation, v singl = 1, g singl = v oct = g oct = 0. They can be calculated by QCD 
sum rules. From our above analysis for the color-singlet operator, we have 

Vsingi - 1 + 0(1(T 2 ), g smgl ~ 0(KT 2 ). (19) 

The sum rules for the color-octet operator are obtained as 

f(a s GG)~ 2/ AT e 4R / f 

f(a s GG)~ 2 . A f e 4A / f (2 ° } 

9 °« = -{-i^ T+m ^ q) )i^n^ B - 

The numerical results for v oct and g oct are given in Fig. 6 from which we obtain 

v oct ~ (1.35 ± 0.76) x 10' 1 , g oct ~ (0.74 ± 0.40) x 10" 1 . (21) 
They are typically at the order of 0.1. 

From the QCD sum rule point of view, the vacuum saturation for the hadronic 
matrix elements of the color-singlet four-quark operators given in Eq. (5) or in the first 
equation of (18) is valid up to a few percent level. However, that for the color-octet 
four-quark operators is violated at ten percent level. The result that the nonfactorizable 
effect for color-octet operators is more significant can be also understood qualitatively, 
if we look at the perturbative diagrams. This effect is at three-loop level, whereas that 
for color-singlet operators at the four-loop level. 

The I /nib corrections to above results can be analyzed in principle. While having 
little influence on our above arguments and calculations, they formally belong to the 
0(1 /ml) effects to the b-hadron lifetimes. 

In Ref . [6] , the extraction of the parameters from B decays are discussed. The role 
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of parameter g S mgi might be significant in comparing the decay rates of B° — > lvX u and 
B~ — > lvX u . From our analysis, the difference of the decay rates should be small due 
to g S ingi is very small even at the scale of m^. The ratio of g S in g i/goct can be obtained 
from the lepton spectrum of the above decays. While having not determined the sign 
of g singh our analysis prefers an even smaller ratio of \g S mgi/ 9oct\ ~ 1/10 than Ref. [6] 
did. These points will be verified by the experiments in the near future. 



The b-hadron lifetime ratios are expressed as [7] 

4f^v = 1 + 0.035! - 0.71 ei + 0.20e 2 + O(^) , 

r\i ' - ( 22 ) 

= 0.98-0.17 ei + 0.20e 2 - (0.013 + 0.0225)r + O(^) . 

In above equation, the energy scale for the parameters is at m b . the B meson decay 
constant fs is taken as 200 MeV. By using the results in Eq. (17), we obtain 

±0.02 and ^ .98 ± 0.01 . (23) 

In obtaining above numbers, the QCD sum rule results for the baryonic parameters r 
and B [8] have been also taken values at scale m&. From Eq. (23) we see that, while 
t(B~)/t(B°) is agree with measurement, the QCD sum rule result for t(A;>)/t(-B ), 
after including 1/mf corrections, still contradicts the experimental measurement. 



In summary, the nonfactorizable contributions of the hadronic matrix elements of 
four-quark operators relevant to the B meson lifetimes have been studied by QCD 
sum rules in the framework of HQET. The vacuum saturation for color-singlet four- 
quark operators is justified at hadronic scale, and the nonfactorizable effect is at a few 
percent level. However, the vacuum saturation for color-octet four-quark operators is 
violated sizably that the nonfactorizable effect cannot be neglected for the B meson 
lifetimes. The implication to the extraction of some of the parameters from B decays 
has been discussed. The B meson lifetime ratio has been predicted to be consistent 
with experiment, t(B~)/t(B°) = 1.09±0.02. However, the discrepancy between theory 

10 



and experiment on t(A{,)/t(_B°) has not been solved. It is unlikely that higher l/m& 
corrections will give the solution. If the experimental data of r(A^) /t(B°) is further 
confirmed in the future, that may imply either the failure of local duality assumption 
or some new physics. 
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Figure captions 



Fig. 1. Dominant non- vanishing Feynman diagrams for T°(u,u'). 

Fig. 2. Nonfactorizable diagrams for a general color-singlet four-quark operator. 

Fig. 3. Nonfactorizable diagrams by two gluon exchanges for T°(uj,uj'). 

Fig. 4. Condensate contribution to T t (uj,uj'). 

Fig. 5. Sum rules for e\ (a) and t<i (b). The sum rule window is T = 0.7— 1.0 GeV. 
Fig. 6. Sum rules for v oct (a) and g oct (b). 
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Fig. 3 
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Fig. 4 
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Fig. 5(a) 




Fig. 5(b) 
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Fig. 6(a) 




Fig. 6(b) 
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